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We have performed a search for new particles which decay to two photons using 1.2 fb~^ of 
integrated luminosity from pp collisions at = 1.96 TeV collected using the CDF II Detector at 



the Fermilab Tevatron. We find the diphoton mass spectrum to be in agreement with the standard 
model expectation, and set limits on the cross section times branching ratio for the Randall-Sundrum 
graviton, as a function of diphoton mass. We subsequently derive lower limits for the graviton mass 
of 230 GeV /c? and 850 GeV/c^, at the 95% confidence level, for coupling parameters (k/Mpt) of 
0.01 and 0.1 respectively. 

PACS numbers: 12.38. Qk 13.85.Rm 13.85.Qk 



A potential signature for new, heavy particles is 
a narrow mass resonance decaying to two energetic 
photons. Such a signature could arise from extra 
spatial dimensions, as in the Randall-Sundrum (RS) 
model [H . In this Letter we present the results of the 
search for this signature in pp collision data from the 
CDF II Detector at the Fermilab Tevatron. 

Some string theories propose that as many as 
seven new spatial dimensions may exist, with their 
geometry potentially responsible for the apparent 
weakness of gravity (the hierarchy problem). In 
the minimal RS model, this geometry is a five- 
dimensional space with negative curvature, bounded 
by two three-dimensional, spatially extended mem- 
branes (or 'branes'). These branes are separated by 
a distance irrc, where Tc is the compactification ra- 
dius. Using (j) as the coordinate of the extra dimen- 
sion (0 < < tt), the standard model (SM) particles 
are confined to the "TeV" brane, located at (f> = it, 
while the gravitational wavefunction is localized at 

= 0. The scale of physical phenomena on the 
TeV brane is specified by A^r = Mpi e^'^^"'^ , where 
Mpi = Mpi/VS^ = 2.4 X IQis GeV is the effective 
four-dimensional (reduced) Planck scale and fc is a 
curvature parameter of the order of the Planck scale. 

To remove the hierarchy between gravity and 
the electroweak force, A^r is set to approximately 

1 TeV. The Planck scale then arises from the small 
overlap of the graviton wave function with the TeV 
brane in the fifth dimension, and the compactifi- 
cation scale gives rise to a Kaluza Klein tower of 
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graviton states, an infinite set of four-dimensional 
particles with increasing masses. This mass spec- 
trum is given by m„ ~ Xn{k/ M pi)Kt^, where Xn are 
the roots of the first-order Bessel function, and the 
states couple with strength 1/A^. Thus, the widths 
and masses of the resonances are dependent on the 
parameter k/Mpi. 

The values of k must be large enough to be con- 
sistent with the apparent weakness of gravity, but 
small enough to prevent the theory from becoming 
non-perturbative Given these considerations, we 
examine values in the range 0.01 < k/Mpi < 0.1. 
For this range, graviton production results in a 
diphoton mass peak narrower than the CDF detector 
resolution. The spin-2 nature of the graviton, decay- 
ing by either s- or p-wave states, favors searches in 
the diphoton channel, where the branching ratio is 
twice that of any single dilepton channel. 

Existing lower mass limits on RS gravitons 
from a search using the D0 detector, in the com- 
bined diphoton, dielectron, and dimuon channels, 
are 250 GeV/c2 for k/Mpi = 0.01 and 785 GeV/c^ 
for k/Mpi = 0.1 at the 95% confidence level 
(C.L.) 3J. The previous limits from the CDF collab- 
oration are from a combined dielectron and dimuon 
search, with limits of 170 GeV/c^ and 710 GeV/c^ 
for fc/Mpi = 0.01 and 0.1 respectively 0]. 

This search uses 1.2 fb~^ of integrated luminos- 
ity collected by the CDF II Detector operating at 
^/s = 1.96 TeV. The detector Q is approximately 
forward-backward and azimuthally symmetric. A 
seven-layer (eight-layer in the forward region) sili- 
con tracker [6| is surrounded by an open-cell drift 
chamber (COT) 0. The fiducial coverage of the 
COT is \t]\ < 1.0, and the silicon detector extends 
the tracking coverage to \ri\ < 2.0. The integrated 
tracking system is contained within a superconduct- 
ing solenoid, providing a 1.4 T magnetic field. Sur- 
rounding these are the electromagnetic and hadronic 
calorimeters [8], divided into "central" (It^I < 1.1) 
and "plug" (1.1 < |?7| < 3.6) regions, providing mea- 
surements of both shower energy and position. At 
the approximate electromagnetic shower maximum, 
the calorimeters contain fine-grained detectors Q 
that measure the shower shape and centroid posi- 
tion in the two dimensions transverse to the shower 



development. Surrounding these detectors is a sys- 
tem of muon detectors [l^l ■ Three levels of real-time 
trigger systems are used to filter events. 

The collision events recorded for analysis were 
selected by at least one of four triggers. Two of 
the triggers require two clusters of electromagnetic 
energy [5|: one of these triggers requires both clus- 
ters to have transverse energy Et > 12 GeV, and 
to be isolated in the calorimeter; while the other re- 
quires the two clusters to have Et > 18 GeV, but 
makes no isolation requirement. To ensure very high 
trigger efficiency at the largest photon Et, events 
are also accepted from two single-photon triggers. 
The first of these triggers requires Et > 50 GeV 
without an isolation criterion; the second requires 
Et > 70 GeV, without an isolation ctriteria, and 
relaxed requirements on the hadronic energy associ- 
ated with cluster. The combination of these triggers 
is effectively 100% efficient for the kinematic region 
used in this search for diphoton events above a mass 
of 50 GeV/c2. 

The triggered events are required to have been 
recorded while the relevant detector elements were 
fully operational. Events are examined in which ei- 
ther both photons are in the central calorimeter (ap- 
proximately in the region \r]\ < 1.04), or one photon 
is in the central and the other is in the plug calorime- 
ter (over the region 1.2 < I77I < 2.8). The data 
were collected between February, 2002 and Febru- 
ary, 2006 and correspond to an integrated luminos- 
ity of 1.2 fb~^ for the central-central (CC) diphoton 
sample and 1.1 fb~^ for the central-plug (CP) sam- 
ple. The luminosity of the CP sample is reduced be- 
cause good silicon detector conditions are required 
for tracking in the plug region. 

In both the CC and CP cases, each of the two 
photons is required to produce a highly electromag- 
netic energy cluster with Et > 15 GeV and to- 
gether to have a reconstructed diphoton invariant 
mass greater than 30 GeV/c^. Both are required 
to be in the fiducial region of the shower maximum 
detectors and to pass the following photon identifi- 
cation requirements: transverse shower profiles con- 
sistent with a single photon, additional transverse 
energy in the calorimeter in a cone of angular ra- 
dius R = y/{A<j)y + (A77)2 = 0.4 around the pho- 
ton candidate < 2 GeV, and the scalar sum of the 
transverse momentum pt of the tracks in the same 
cone < 2 GeV/c. Photons in the central detector are 
required to produce isolated energy clusters in the 
shower maximum detector. 

The selected data consist of 11, 088 CC and 
20, 933 CP events. The diphoton invariant mass 
distribution for these events, histogrammed in bins 



equivalent to the mass resolution (approximated by 
0.13^m(GcV/c2) e 0.014m GeV, where m is dipho- 
ton mass) is shown in Fig. [TJ The highest mass 
pairs occur at 602 GeV/c^ for the CC data and 
454 GeV/c^ for the CP data, as shown in Fig. H 

The expected numbers of RS graviton events, as 
a function of graviton mass, are estimated using the 
HERWIG 6.510 event generator [13], with CTEQ5l 
parton distribution functions (PDFs) [l^], and pro- 
cessed by the GEANT 3 based CDF II detector simu- 
lation [lj|. The combined acceptance and selection 
efficiency for RS events increases from 0.31±0.01stat 
for gravitons of mass 200 GeV/c^ to 0.36 ± O.Olsfat 
for gravitons of mass 1050 GeV/c^. Final corrections 
to the efficiency are derived by comparing the mea- 
sured and Monte Carlo simulated detector response 
to electrons from Z boson decays, since a pure sam- 
ple of reconstructed photons is not available, and the 
characteristics of energy deposited in the calorimeter 
by electrons are almost indistinguishable from those 
of photons. The largest systematic uncertainties on 
the expected number of gravitons arise from the lu- 
minosity measurement (6%) and the choice of PDF 
(4%). The latter uncertainty is determined from the 
variation in the efficiency when employing different 
PDF parameterizations. 

Most Z bosons that decay to e+e" are rejected 
by requiring no associated tracks; however, approx- 
imately 1% of these remain in the CP sample at a 
mass peak below the search region. There are two 
other significant components to the diphoton data 
sample. The first is SM diphoton production, which 
is estimated with the diphox ncxt-to-leading-order 
(NLO) Monte Carlo [13] calculation. This program 
computes a cross section as a function of mass. The 
mass spectrum is then corrected by an efficiency 
function derived from a SM diphoton sample gener- 
ated by PYTHIA [l6l| and processed through the full 
detector simulation. The leading systematic uncer- 
tainty on this background estimate is approximately 
20%, which comes from the variation in the cross 
section when the renormalization scale, Q^, is var- 
ied in the generation. The second component is jets, 
consisting of quarks that fragment into a high mo- 
mentum tt", which subsequently decays as tt^ — > 77. 
The resulting photon showers may overlap, and can 
pass the photon selection. The diphoton mass distri- 
bution of this background is derived from a sample 
of photon-like jets obtained by loosening the photon 
selection criteria, then removing the events which 
pass all the signal selection requirements. The selec- 
tion of the photon-like jets is varied and the result- 
ing variation in the shape of the mass distribution is 
approximately 20%. This is the leading systematic 



uncertainty on the jet background. 

To extrapolate to the highest masses, the mass 
distributions of the two background sources are fit- 
ted to a product of a polynomial and the sum of two 
(for photon-like jets) or three (for diphox) exponen- 
tials. The fit to the photon-like jets includes a con- 
tribution from the DIPHOX shape to allow for true 
diphoton events which appear in this sample. For 
the signal region background estimate, the DiPHOX 
distribution is normalized to the NLO cross section 
predicted by diphox. The jet background shape is 
normalized to the observed shape in the invariant 
mass region 30 GeV/c^ to 100 GeV/c^ after sub- 
tracting the estimated background from SM dipho- 
ton production. 

Figure [2] shows the observed mass spectra com- 
pared to the predictions. This estimate of the back- 
ground, which shows good agreement with the obser- 
vation, is only for an a priori comparison, and is not 
used in setting upper limits on a specific production 
model. 

To find the most accurate description of the 
background for setting limits, we fit the CC and 
CP mass spectra to a polynomial multiplied by a 
sum of three components: two exponentials and the 
DIPHOX shape. The normalization of each compo- 
nent is allowed to float in the fit. This function is 
used for the background estimate as a function of 
mass with statistical uncertainties propagated from 
the fits. DIPHOX systematic uncertainties, as de- 
scribed above, are also propagated. 

The background shape is combined with the 
simulated signal shape in a binned likelihood 
method, with the contents of the bins treated with 
Poisson statistics. The likelihood as a function of 
cross section is given by 



L{a) = n 



-Mi('T) 



where 



N 



(1) 



(2) 



tot 



Here Nf , and Nf are respectively the num- 
bers of observed, background and signal events in 
bin I, and N^^^ is the total number of signal events 
passing selection requirements. The parameter e is 
the total acceptance and efficiency for selecting the 
events predicted by the model, C is the integrated 
luminosity, and a is the model cross section. Sys- 
tematic uncertainties are included in the likelihood 
with a gaussian Bayesian prior. 



k/Mp, 


Lower Mass Limit (GeV/c^) 


0.1 


850 


0.07 


784 


0.05 


694 


0.025 


500 


0.01 


230 



TABLE I: The 95% C.L. lower limits on the mass of the 
RS graviton for the specified values of k/Mpi. 



The likelihood function is calculated separately 
for the CC and CP channels, then a combined likeli- 
hood is formed from the product of the two, account- 
ing for the correlations among systematic uncertain- 
ties. We obtain a posterior density in a assuming a 
uniform prior in a applying Bayes' Theorem. The 
95% C.L. upper limit corresponds to the point be- 
low which 95% of the integral of the posterior density 
lies. 

The result is shown in Fig. [31 as a function 
of graviton mass, along with the theoretical cross 
section times branching ratio for RS gravitons with 
k/Hpi set to 0.1, 0.07, 0.05, 0.025 and 0.01. The 
leading-order graviton production cross section is 
multiplied by a factor of 1.3, to correct for dia- 
grams at higher-order in as |17l |. From the limit 
on a ■ B{G 77), lower mass bounds are obtained 
for the first excited state of the RS graviton as a 
function of the parameter k/Mpi. The 95% C.L. 
excluded region in the k/Mpi and graviton mass 
plane is displayed in Fig. 31 with the mass limits 
summarized in Table [H 



In conclusion, we have searched for evidence 
of an anomalous peak in the diphoton mass spec- 
trum using approximately 1.2 fb~^ of data collected 
by the CDF II Detector at the Fermilab Tevatron. 
Our sample includes events with both photons in 
the central detector, and events with one photon in 
the central detector and one in the plug detector. 
We find no evidence of new physics. We evaluate 
one model of hypothetical new diphoton production 
and exclude RS gravitons below masses ranging from 
230 GeV/c^ to 850 GeV/c^, for a coupling parameter 
k/Hpi of 0.01 to 0.1, at the 95% C.L. This results 
in a significant improvement, at high mass, over the 
previous best available limit. 
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FIG. 1: The diphoton invariant mass distribution of 
events for both CC and CP channels, histogrammed in 
bins of approximately one unit of calorimeter mass res- 
olution. The enhancement near 90 GeV /(? is due to 
misidentified Z boson decays to electrons in the CP sam- 
ple. 




Total background 
Jets faking photons 



a) -i 



400 500 

m(yY) (GeV/c ) 



600 700 




FIG. 2: The mass distribution in the CC (a) and CP 
(b) signal regions with the a priori background overlaid. 
The points are the data. The dotted line shows the jets 
which fake photons as predicted from the photon-like jet 
sample, and the solid line shows this background plus 
the DIPHOX SM diphoton distribution. The grey band 
shows the uncertainty on the total background. This 
background is not used in setting limits. 




FIG. 3: The 95% C.L. upper limit on the production 
cross section multiplied by branching fraction of an RS 
model graviton decaying to diphotons (it • B{G —> 77)), 
as a function of graviton mass. Also shown are the pre- 
dicted (cr • B) curves for fc/Mp; = 0.01, 0.07, 0.05, 0.025 
and 0.1. 




FIG. 4: The 95% C.L. excluded region in the plane of 
k/Mpi and graviton mass. 



